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We place direct upper limits on the amplitude of gravitational waves from 28 isolated radio 
pulsars by a coherent multi-detector analysis of the data collected during the second science run of 
the LIGO interferometric detectors. These are the first direct upper limits for 26 of the 28 pulsars. 
We use coordinated radio observations for the first time to build radio-guided phase templates for 
the expected gravitational wave signals. The unprecedented sensitivity of the detectors allow us to 
set strain upper limits as low as a few times 10~ 24 . These strain limits translate into limits on the 
equatorial ellipticities of the pulsars, which are smaller than 10 for the four closest pulsars. 
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A worldwide effort is underway to detect gravitational 
waves (GWs) and thus test a fundamental prediction of 
General Relativity. In preparation for long-term opera- 
tions, the LIGO and GEO experiments conducted their 
first science run (SI) during 17 days in 2002. The de- 
tectors and the analyses of the SI data are described in 
Refs. 0] and 0-Q , respectively. LIGO's second science 
run (S2) was carried out from 14 Feb to 14 April 2003, 
with dramatically improved sensitivity compared to SI. 
During S2 the GEO detector was not operating. 

A spinning neutron star is expected to emit GWs if it 
is not perfectly symmetric about its rotation axis. The 
strain amplitude ho of the emitted signal is proportional 
to the star's deformation as measured by its ellipticity e 
Using data from S2, this paper reports direct obser- 
vational limits on the GW emission and corresponding 
cllipticities from the 28 most rapidly rotating isolated 
pulsars for which radio data is complete enough to guide 
the phase of our filters with sufficient precision. These 
are the first such limits for 26 of the pulsars. We con- 
centrate on isolated pulsars with known phase evolutions 
and sky positions to ensure that our targeted search re- 
quires relatively few unknown parameters. 

The limits reported here are still well above the indi- 
rect limits inferred from observed pulsar spindown, where 
available (Fig.^). However, fourteen of our pulsars are in 
globular clusters, where local gravitational accelerations 
produce Dopplcr effects that mask the intrinsic pulsar 
spindown, sometimes even producing apparent spinup. 
For these pulsars our observations therefore place the first 
limits that are inherently independent of cluster dynam- 
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ics, albeit at levels well above what one would expect if 
all globular cluster pulsars are similar to field pulsars. 

Our most stringent ellipticity upper limit is 4.5 x 10 -6 . 
While still above the maximum expected from conven- 
tional models of nuclear matter, distortions of this size 
would be permitted within at least one exotic theory of 
neutron star structure 0. 

Detectors. — LIGO comprises three detectors. Each de- 
tector is a power-recycled Michelson interferometer, with 
Fabry-Perot cavities in the long arms. A passing GW 
produces a time- varying differential strain in these arms, 
and the resulting differential displacement of the cavity 
test mass mirrors is sensed intcrfcromctrically. Two de- 
tectors, the 4 km-arm HI and the 2 km-arm H2 detectors, 
are collocated in Hanford WA. The 4 km-arm LI detec- 
tor is situated in Livingston Parish LA. Improvements 
in noise performance between SI and S2 were approx- 
imately an order of magnitude over a broad frequency 
range. Modifications that were made between SI and S2 
to aid in noise reduction and improve stability include: 

i) increased laser power to reduce high-frequency noise, 

ii) better angular control of the mirrors of the interfer- 
ometer and iii) the use of lower noise digital test mass 
suspension controllers in all detectors. 

During S2, the LIGO detectors' noise performance in 
the band 40-2000 Hz was better than any previous inter- 
ferometer. The best strain sensitivity, achieved by LI, 
was ~3x 10 -22 Hz -1 ^ 2 near 200 Hz (where it translates 
through Eq. (2.2) of into a detectable amplitude for a 
continuous signal of about 3 x 10~ 24 , as shown in Fig.QJ. 
The relative timing stability between the interferometers 
was also significantly improved. Monitored with GPS- 
synchronized clocks to be better than 10 /is over S2, it 
allowed the coherent combination of the strain data of 
all three detectors to forrnjoint upper limits. 

Analysis method. — In [2] a search for gravitational 
waves from the millisecond pulsar J1939+2134 using SI 
data was presented. In that work, two different data 
analysis methods were used, one in the time domain and 
the other in the frequency domain. Here we extend the 
former method 00 and apply it to 28 isolated pulsars. 

Following we model the sources as non-precessing 
triaxial neutron stars showing the same rotational phase 
evolution as is present in the radio signal and perform a 
complex heterodyne of the strain data from each detec- 
tor at the instantaneous frequency of the expected grav- 
itational wave signal, which is twice the observed radio 
rotation frequency. These data are then down-sampled 
to 1/60 Hz and will be referred to as Bf.. Any gravi- 
tational signal in the data would show a residual time 
evolution reflecting the antenna pattern of the detector, 
varying over the day as the source moved through the 
pattern, but with a functional form that depended on 
several other source-observer parameters: the antenna 
responses to plus and cross polarisations, the amplitude 
of the gravitational wave ho, the angle between the line- 
of-sight to the pulsar and its spin axis t, the polarisation 
angle of the gravitational radiation ijj (all defined in |6j) 
and the phase <po of the gravitational wave signal at some 
fiducial time to- Let a be a vector in parameter space 



3 




x 10 



10' 10" 
Frequency (Hz) 

FIG. 1: Upper curves: ho amplitudes detectable from a known 
generic source with a 1% false alarm rate and 10% false dis- 
missal rate, as given by Eq. (2.2) in for single detector 
analyses and for a joint detector analysis. All the curves 
use typical S2 sensitivities and observation times. HI and 
H2 are the 4km-arm and the 2km-arm detectors located in 
Hanford WA. LI is the 4km-arm detector situated in Liv- 
ingston Parish LA. Lower curve: LIGO design sensitivity for 
1 yr of data. Stars: upper limits found in this paper for 28 
known pulsars. Circles: spindown upper limits for the pulsars 
with negative frequency derivative values if all the measured 
rotational energy loss were due to gravitational waves and 
assuming a moment of inertia of 10 45 g cm 2 . 



with components (ho, t, ip, 4>a)- 

The analysis proceeds by determining the posterior 
probability distribution function (pdf) of a given the 
data Bk and the signal model: 



p(a\{B k }) cx p(a)p({B k }\a) 



(1) 



where p({B k }\a) is the likelihood and p(a) the prior pdf 
we assign to the model parameters. We have used a uni- 
form prior for cos/, <Po. and ho (ho > 0), in common 
with [2| • A uniform prior for ho has been chosen for its 
simplicity and so that our results can readily be compared 
with other observations. This prior favors high values of 
ho (which comprise the majority of the parameter space) 
and therefore generates a somewhat conservative upper 
limit for its value. Indeed the reader might prefer to re- 
gard our resulting posterior pdfs for ho as marginalised 
likelihoods rather than probabilities for ho — these are 
functionally identical using our priors. 

As in Q we use a Gaussian joint likelihood for 
p({Bk}\a). In the SI noise floor was estimated over 
a 60s period from a 4 Hz band about the expected sig- 
nal frequency. This gave a reliable point estimate for the 
noise level but was sensitive to spectral contamination 
within the band, as demonstrated in the analysis of the 
GEO SI data. In this paper we exploit the improved sta- 
tionarity of the instruments that make it reasonable to 
assume the noise floor is constant over periods of 30 min 
In addition we restrict the bandwidth to 1/60 Hz, 
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FIG. 2: Parameters of the artificial pulsar PI, recovered from 
12 h of strain data from the Hanford and Livingston inter- 
ferometers. The results are displayed as marginal pdfs for 
each of the four signal parameters. The vertical dotted lines 
show the values used to generate the signal, the colored lines 
show the results from the individual detectors (HI green, H2 
blue, LI red), and the black lines show the joint result from 
combining coherently data from all three. 



which makes it possible to search for signals from pul- 
sars at frequencies close to strong spectral disturbances. 
However, the noise level now determined is less certain 
as the estimate relies on fewer data. We take account 
of this increased uncertainty by explicitly marginalising 
with a Jeffreys prior over the constant but unknown noise 
level for each 30 min period of data . The likelihood for 
this analysis is then the combined likelihood for all the 
30 min stretches of data, labeled by j in Eq. , taken 
as independent: 

p({B k }\a) cx UpdB^la), (2) 

3 

I k 2{f) 

p({B k } 3 \a) (x J2 \ B k~Vk\ 2 | 



fe=fci 



where y& is the signal model given by Eq. (4.10) in @ 



and m = k 



2(i) 



fci 



1 = 30 is the number of B k data 



points in a 30 min segment. 

In principle the period over which the data are as- 
sumed stationary need not be fixed, and can be adjusted 
dynamically to reflect instrumental performance over the 
run. We have limited our analysis to continuous 30 min 
stretches of data, which included more than 88% of the 
S2 science data set. Inclusion of shorter data sections 
would at best have resulted in a ~ 6% improvement on 
the strain upper limits reported here. 

Validation by hardware injections. — All the soft- 
ware used in this analysis is available in the 
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LAL and LALapps CVS repositories (www.lsc- 
group.phys.uwm.edu/daswg/projects/lal.html) with 
the tag "pulgroup_paper_s2tds" . It was validated by 
checking its performance on fake pulsar signals injected 
in artificial and real detector noise, both in software (Q) 
and in hardware. In particular, two artificial signals (PI, 
P2) were injected into all three detectors by modulating 
the mirror positions via the actuation control signals 
with the strain signal we should expect from a hypo- 
thetical pulsar. These injections were designed to give 
an end-to-end validation of the search pipeline starting 
from as far up the observing chain as possible. 

The pulsar signals were injected for 12 h at frequen- 
cies of 1279.123 Hz (PI) and 1288.901Hz (P2) with fre- 
quency derivatives of zero and — 10 -8 Hz s" 1 respectively, 
and strain amplitudes of 2 x 10~ 21 . This gives signal-to- 
noise ratios (as defined by Eq. (79) of |||) of 26 and 40 
for PI in HI and LI respectively and of 38 and 34 for 
P2. The signals were modulated and Doppler shifted to 
simulate sources at fixed positions on the sky with "0 = 0, 
cos l = and <fio — 0. To illustrate, posterior pdfs for the 
recovered PI signal are shown in Fig. [2 The results de- 
rived from the different detectors are in broad statistical 
agreement, confirming that the relative calibrations are 
consistent and that the assessments of uncertainty (ex- 
pressed in the posterior widths) are reasonable. Results 
for P2 were very similar to these. 

The phase stability of the detectors in S2 allowed us to 
implement a joint coherent analysis based on data from 
all three participating instruments. This technique was 
noted in @], but could not be performed on the SI data 
because of timing uncertainties that existed when those 
observations were performed. The black lines in Fig. [2] 
show marginalisations of the joint posterior from HI, H2 
and LI, i.e., 

p(a|Hl,H2,Ll) oc jj(a)jj(Hl|a)jj(H2|o)p(U|a). (4) 

In an ideal case of three detectors of similar sensitivi- 
ties and operational times these coherent results would 
be approximately times tighter than the individual 
results. The posteriors for <j>Q clearly highlight the rela- 
tive coherence between the instruments and verify that 
similar joint methods can be used to set upper limits on 
our target pulsars. 

Results. — From the ATNF pulsar catalogue 
(www. atnf . csiro . au/research/pulsar/psrcat/) 
we selected 28 isolated pulsars with rotational frequen- 
cies greater than 20 Hz and for which sufficiently good 
timing data were available (Table I). For 18 of these, we 
obtained updated timing solutions from regular timing 
observations made at the Jodrell Bank Observatory 
using the Lovell and the Parkes telescopes, adjusted for 
a reference epoch centred on the epoch of the S2 run 
(starred pulsars in Table I). Details of the techniques 
that were used to do this can be found in [ljj. We also 
checked that none of these pulsars exhibited a glitch 
during this period. 

The list includes globular cluster pulsars (including 
isolated pulsars in 47 Tuc and NGC6752), the SI tar- 
get millisecond pulsar (J1939+2134) and the Crab pul- 
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1(1 


2.4 


1.8 


J1024-0719* 


193.72 


-6.95x10" 


16 


3.9 


0.86 


B1516+02A 


180.06 


-1.34x10" 


15 


3.6 


21 


J1629-6902 


166.65 


-2.78x10" 


16 


2.3 


2.7 


J1721-2457 


285.99 


-4.80x10" 


16 


4.0 


1.8 


J1730-2304* 


123.11 


-3.06x10" 


16 


3.1 


2.5 


J1744-1134* 


245.43 


-5.40x10" 


16 


5.9 


0.83 


J1748-2446C 


118.54 


+8.52x10" 


15 


3.1 


24 


B1820-30A* 


183.82 


-1.14x10" 


13 


4.2 


24 


B1821-24* 


327.41 


-1.74x10" 


13 


5.6 


7.1 


J1910-5959B 


119.65 


+1.14x10" 


14 


2.4 


8.5 


J1910-5959C 


189.49 


-7.90x10" 


17 


3.3 


4.7 


J1910-5959D 


110.68 


-1.18x10" 


14 


1.7 


7.2 


J1910-5959E 


218.73 


+2.09x10" 


14 


7.5 


7.9 


J1913+1011* 


27.85 


-2.61x10" 


12 


51 


6 900 


J1939+2134* 


641.93 


-4.33x10" 


14 


13 


2.7 


B1951+32* 


25.30 


-3.74x10" 


12 


48 


4 400 


J2124-3358* 


202.79 


-8.45x10" 


16 


3.1 


0.45 


J2322+2057* 


207.97 


-4.20x10" 


16 


4.1 


1.8 



TABLE I: The 28 pulsars targeted in the S2 run, with approx- 
imate spin parameters. Pulsars for which radio timing data 
were taken over the S2 period are starred (*). The right- 
hand two columns show the 95% upper limit on ho, based on 
a coherent analysis using all the S2 data, and corresponding 
ellipticity values (e, see text). These upper limit values do 
not include the uncertainties due to calibration and to pulsar 
timing accuracy, which are discussed in the text, nor uncer- 
tainties in the pulsar's distance, r. 



sar (B0531+21). Although Table I only shows approx- 
imate pulsar frequencies and frequency derivatives, fur- 
ther phase corrections were made for pulsars with mea- 
sured second derivatives of frequency. Timing solutions 
for the Crab were taken from the Jodrell Bank online 
ephemeris [Til ], and adjustments were made to its phase 
over the period of S2 using the method of [l2^ . 

The analysis used 910 hours of data from HI, 691 hours 
from H2, and 342 hours from LI. There was no evidence 
of strong spectral contamination in any of the bands in- 
vestigated, such as might be caused by an instrumen- 
tal feature or a potentially detectable pulsar signal. A 
strong gravitational signal would generate a parameter 
pdf prominently peaked off zero with respect to its width, 
as for the hardware injections. Such a pdf would trigger 
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a more detailed investigation of the pulsar in question. 
No such triggers occurred in the analysis of these data, 
and we therefore present upper limits. 

The upper limits are presented as the value of ho 
bounding 95% of the cumulative probability of the 
marginalised strain pdf from ho = 0. The joint upper 



limit ft,Q 5% therefore satisfies 



0.95 = 



dh 



hn=0 



p(a|Hl,H2,Ll)dtdV>d<fo, (5) 



consistent with The uncertainty in the noise floor 
estimate is already included, as outlined above. 

The remaining uncertainties in the upper limit values 
of Table I stem from the calibration of the instrument 
and from the accuracy of the pulsar timing models. For 
LI and H2, the amplitude calibration uncertainties are 
conservatively estimated to be 10% and 8%, respectively. 
For HI, the maximum calibration uncertainty is 18%, 
with typical values at the 6% level. Phase calibration 
uncertainties are negligible in comparison: less than 10° 
in all detectors. Biases due to pulsar timing errors are 
estimated to be 3% or less for J0030+0451, and 1% or 
less for the remaining pulsars (see 0] for a discussion of 
the effect of these uncertainties). 

Discussion. — The improved sensitivity of the LIGO 
interferometers is clear from the strain upper limit for 
PSR J1939+2134, which is more than a factor of ten 
lower than was achieved with the SI data Q ■ In this anal- 
ysis the lowest limit is achieved for PSR J1910— 5959D at 
the level of 1.7 x 10 -24 , largely reflecting the lower noise 
floor around 200 Hz. 

Table I also gives approximate limits to the cllipticitics 
of these pulsars from the simple quadrupole model 



0.237 



ho 



r 1 Hz 2 10 45 g cm 2 



10- 24 lkpc f 2 



(6) 



where r is the pulsar's distance, which we take as the 
dispersion measure distance using the model of Taylor 



and Cordes |T3 |. and I zz its principal moment of inertia 
about the rotation axis, which we take as 10 45 gcm 2 . 

As expected, none of these upper limits improves on 
those inferred from simple arguments based on the gravi- 
tational luminosities achievable from the observed loss of 
pulsar rotational kinetic energy. However, as discussed 
in the introduction, for pulsars in globular clusters such 
arguments are complicated by cluster dynamics, which 
the direct limits presented here avoid. 

The result for the Crab pulsar (B0531+21) is within 
a factor of about 30 of the spindown limit and over an 
order of ma gnit ude better than the previous direct up- 
per limit of |l4j . The equatorial ellipticities of the four 
closest pulsars (J0030+0451, J2124+3358, J1024-0719, 
and J1744— 1134) are constrained to be less than 10 -5 . 

Once the detectors operate at design sensitivity for a 
year, the observational upper limits will improve by more 
than an order of magnitude. The present analysis will 
also be extended to include pulsars in binary systems, 
significantly increasing the population of objects under 
inspection. 
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